Cooperativity and Frustration in Protein-Mediated Parallel Actin Bundles 
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We examine the mechanism of bundling of cytoskeletal actin filaments by two representative 
bundling proteins, fascin and espin. Small-angle X-ray studies show that increased binding from 
linkers drives a systematic overtwist of actin filaments from their native state, which occurs in a 
linker-dependent fashion. Fascin bundles actin into a continuous spectrum of intermediate twist 
states, while espin only allows for untwisted actin filaments and fully-overtwisted bundles. Based on 
a coarse-grained, statistical model of protein binding, we show that the interplay between binding 
geometry and the intrinsic flexibility of linkers mediates cooperative binding in the bundle. We 
attribute the respective continuous/discontinous bundling mechanisms of fascin/espin to differences 
in the stiffness of linker bonds themselves. 

PACS numbers: 



Actin binding proteins (ABP) that direct the assem- 
bly of F-actin cytoskeletal polymers are often divided 
into two classes, those that induce formation of net- 
works, and those that induce formation of finite-sized 
parallel bundles [lj. These motifs have been observed 
for a variety of linkers, from ABP's to simple multivalent 
ions, and have been studied theoretically and experimen- 
tally @, i |, S, @, 0, & & M EL El- Espin and fascin 
are two representative bundle-forming ABP's. Espins 
are found in mechanosensory microvilli and microvillar 
derivatives, while fascin is typically found in filopodia. 
Although the gross structure of the induced F-actin bun- 
dles are similar for espin and fascin 0, 0| , they behave 
differently, and serve cellular functions with different re- 
quirements. Here, we aim to explore a deeper taxon- 
omy governing the different behaviors of bundle-forming 
ABP's. 

In this Letter we demonstrate that while different 
crosslinkers ultimately drive parallel actin bundles to the 
same structural state, the thermodynamic transition to 
that state depends sensitively on linker stiffness. Mon- 
itoring the structural evolution of bundled filaments by 
Small- Angle X-ray Scattering (SAXS), we find that in- 
creasing the ratio of fascin to actin leads to a continuous 
overtwisting of filaments from their native symmetry In 
contrast, crosslinking by espin produces a coexistence of 
two populations, one with the fully overtwisted geometry, 
and one with native twist. We propose a coarse-grained 
lattice model of crosslinking in actin bundles to capture 
the interplay between filament and crosslinker flexibil- 
ity as well the incommensurate geometries of actin fila- 
ments and fully crosslinked bundles. This model reveals 
that stiffness of crosslinking bonds and resistance to fil- 
ament torsion sensitively control the level of cooperativ- 
ity of crosslinking at different points along the filament. 



The mean-field thermodynamics of this model predicts: 
1) a flexible linker regime allows a continuous increase 
of crosslinks with increased chemical potential; 2) a stiff 
linker regime exhibits a highly cooperative and discontin- 
uous linker binding transition; and 3) a critical-end point 
separating these regimes. The respective continuous and 
discontinuous changes in filament overtwist measured by 
scattering can be correlated with the flexible linker and 
stiff linker regimes of the lattice model, where a similar 
response to increased crosslinking is predicted, suggesting 
that a small differences in linker structure lead qualitative 
differences in global phase behavior of the cytoskeleton. 

To prepare X-ray samples, fresh F-actin and crosslink- 
ing protein were mixed at specific molar ratios R = 



Ncrossiinker/ Nq -actin, with 0.15 mg F-actin. F-actin was 
prepared from rabbit skeletal muscle G-actin monomer 
(Cytoskeleton, Inc.) which was first resuspended in 5 
mM Tris, 0.2 mM CaCl 2 , 0.5 mM ATP, 0.2 mM dithotrc- 
ithol and 0.01% NaN3, pH 8.0 and then polymerized into 
F-actin by adding 100 mM KC1. F-actin was then treated 
with human plasma gelsolin (Cytoskeleton, Inc.) to con- 
trol average F-actin length (~ lum) and with phalloidin 
to prevent depolymerization [13J. The F-actin solution 
was then centrifuged at 100 000 xg for 1 hour to remove 
polymerization buffer and resuspended in E-buffer: 0.1 
M KC1, 10 mM HEPES, 1 mM dithotreithol, 1.5 mM 
NaN3, pH 7.4. Crosslinking proteins included recombi- 
nant rat espin 3A (34.3 kDa) and recombinant human 
fascin (57.8 kDa), which were expressed in bacteria with 
an N-tcrminal 6xHis tag, affinity purified under non- 
denaturing conditions and dialyzed into E-buffer. Sam- 
ples of F-actin mixed with crosslinker were mixed, incu- 
bated, and centrifuged in sealed quartz capillaries. SAXS 
experiments were performed at 9 KeV at beam line 4-2 
of the Stanford Synchrotron Radiation Lightsource and 
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FIG. 1: Experimental evidence of first and second order twisting transitions. 2D SAXS images of (a) fascin-actin bundles and 
(b) espin-actin bundles. Circularly averaged SAXS data showing first and second layer line peaks for (c) fascin-actin bundles 
and (d) espin-actin bundlesQ as a function of R. Data is shown with a Pseudo-Voigt background subtraction. Arrows show 
position of first layer line peak maximum in (c) and the position of the unbundled layerline peaks at 0.114 and 0.125A -1 . (e) 
Measured twist of actin bundles as a function of R. 



at 12 KcV at the BESSRC-CAT (beam line 12-ID) at the 
Advanced Photon Source. The scattered radiation was 
collected using an MAR Research CCD camera (pixel 
size = 79 /im). The sample-to-detector distances are 
set such that the q-range is 0.01 < q < 0.2 A -1 , where 
q = (An sm8)/ \, A is the x-ray wavelength, and 29 is the 
scattering angle. The 2D SAXS data from both beam- 
lines have been checked for mutual consistency. As de- 
scribed previously Q, the twist of the actin filaments 
when bundled with crosslinking proteins was determined 
by fitting 2D SAXS data to the four sphere model of 
variably-twisted F-actin convolved with the bundle struc- 
ture factor [3, 0] • 

The structure of espin-actin and fascin-actin bundles 
has been previously investigated @, 0, EI EE E3, al- 
though the thermodynamic phase behavior of these actin 
+ABP systems has not been mapped out. SAXS data 
for F-actin condensed by fascin or cspin are presented 
in Fig. [TJ The circularly averaged peak positions of 
the hcxagonally coordinated fascin-actin bundle are sim- 
ilar to those of the espin-actin bundles. Peaks were 
found at 0.057, 0.100, 0.120, 0.134 A" 1 for both es- 
pin and fascin mediated bundles at high R, with the 
first two corresponding to the inter-actin structure fac- 
tor peaks, and the latter corresponding to intra- actin he- 
lical layer line peaks. The inter-actin spacing for the 
fascin-actin bundles obtained from the position of the 
qio peak (most intense peak visible in Fig. [T£i,b, was 
equal to 4tt/(v / 39io) = 12.9 ±0.3 nm, slightly larger than 
that of espin-actin bundles, 12.6 ± 0.2 nm [7J. This cor- 
responds to a fascin size of 5.4 ± 0.3 nm, and an espin 
size of 5.1 ± 0.2 nm using an F-actin diameter of 7.5 
nm [3]. Using the 4-spherc model, we found that the 
position of the espin-actin bundle layer lines indicated a 
F-actin overtwist of 0.9 ±0.2 degrees from the native left- 
handed 13/6 monomers/turn twist symmetry of unbun- 
dled F-actin to a symmetry of 28/13 monomers/turn 
Furthermore, at low R, coexisting bundled and unbun- 
dled phases are observed in the espin-actin system, as in 
Fig. [TJi at R = 0.05 where broad 13/6 layer line peaks at 



0.114 and 0.125 A^ 1 can be observed simultaneously with 
sharp Gaussian peaks of the overtwisted hexagonal bun- 
dles. This 2-phase coexistence in the espin-actin bundle 
data, and constant layer line peak position is in strong 
contrast to the small, systematic shift of the first layer 
line peak of the fascin-actin bundles towards higher q ob- 
served with increasing R (Fig. [T]c). This systematic shift 
in peak position is only visible in the layer line peaks of 
the fascin-actin system, and not in the inter-actin struc- 
ture factor peaks, indicating that it is the F-actin twist 
which is gradually increasing from the native F-actin un- 
bundled twist symmetry (13/6 monomers/ turn) with in- 
creasing fascin concentration, with a maximum of ~ 0.9 
degrees of over-twist at high R, in agreement with recent 
measurements @, E3 ■ The contrast in twisting behavior 
for espin and fascin mediated bundles is summarized in 
Fig. QJ. This fascin-actin bundle data shows a similar 
decrease in twist with decreasing fascin concentration to 
that previously published Q . Espin-crosslinked actin ex- 
hibits a jump between coexisting "low" and "high" twist 
states with increasing espin concentration via a first order 
transition, while fascin-crosslinkcd actin exhibits gradual 
twist changes from "low" to "high" twist with increasing 
fascin concentration, suggesting a continuous thermody- 
namic pathway. 

To study cross-linking thermodynamics, we introduce 
a coarse-grained lattice model reflecting the intrinsic geo- 
metrical frustration of parallel actin bundles. The helical 
axes of actin filaments are positioned at the vertices of 
a hexagonal lattice with spacing D. The helical con- 
figuration of G-actin monomers in each filament is de- 
scribed by a set of XY "spins" on planes spaced along 
the backbone of the filament, as pictured in Fig.O The 
ith monomer is then represented by a spin vector aS,, 
where a is the monomer size. The positions described by 
these vectors serve as a proxy for the locations of bind- 
ing sites on the monomers themselves. For the native 
configuration of actin filaments with 13/6 symmetry, the 
spins wind around the filament axis by ujq = 127r/13 per 
monomer (see Fig. HJa)). 
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We introduce a Hamiltonian, described by a set of spins 
for filament configurations, as well as the binding site 
occupancy between two monomers, i and j, separated 
by a distance A i3 -, 



TLb 



ending 



= J2n ij [-e + -(A ij -A ) 



(1) 



where the sum runs over sites on neighboring filaments 
and riij equals or 1. Here eo describes the minimum 
(distortion free) energy of optimal binding and k the stiff- 
ness of linkers, defining the energy cost to deform the 
linkers from an aligned state, Ay = Ao. In this model, 
we consider the in-plane crosslinks so that the square- 
deformation has a rather simple form, (A^ — Ao) 2 ~ 
Co - 2a 2 D.y ■ (Si - Sj) + C(S 2 ), where Dj, is the unit vec- 
tor of a lattice direction. Hence, protein crosslinking oc- 
curs more favorably where two monomers co-orient along 
the directions of the lattice packing. Based on this model 
we predict that a unique regular structure maximizes the 
number of perfectly aligned crosslinks/monomer in the 
bundle, while requiring minimal distortion of the intrin- 
sic twist of the filament [2(| . The structure is composed 
of 4 sections of 5-monomcrs with 30/14 (overtwisted) 
symmetry and 2 sections of 4- monomers with 24/11 (un- 
dertwistcd) symmetry, so that 6 monomers/repeat align 
perfectly with six-fold lattice directions and all bonds 
from neighboring filaments are coincident. This struc- 
ture fulfills an overall repeat unit of 28 monomers per 13 
turns, consistent with the overtwisted 28/13 geometry 
observed by scattering. Based on an exhaustive search, 
we have found that hexagonal tilings of alternative com- 
posite structures of up 40 monomers/repeat - including 
the corresponding 13/6 structure have smaller fraction 
of bound monomers than 6/28 ~ 0.214, provided by the 
composite 28/13 geometry 

The conformational adjustments of filaments required 
for optimal binding give rise to cooperative crosslinking, 
mediated by torsional fluctuations of filament and linker 
flexibility. To demonstrate this, we adopt a continuum 
model for twist distortions, given by angular deviations 
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FIG. 2: (a) A schematic picture of two actin filaments linked 
by a crosslink, (b) G-actin monomers in filaments are repre- 
sented by a set of XY spins, (c) The top view of the angular 
distribution of crosslinkers (red arrows) in the unit cell of 
28/13 groundstate. 



from the native filament geometry, TIq — ^ Y^ti^-fit ~ 
wo) 2 , where £ denotes the vertical layer, C is the tor- 
sional stiffness and A(j)g = <f>£ — §t-\ is the azimuthal 
angle difference between two adjacent monomers along 
the filament. Based on the geometric distortion of bonds 
the energy for adding a bond at a layer I can be writ- 
ten as, — e + U(l — cos[<^ — <f> m ]). Here, (£, m) label the 
vertical and angular position of bonds, and 4> m = 2-Km/Q 
indicates the preferred 6-fold direction of monomer ori- 
entation. U is a measure of linker flexibility, U ~ ka 2 . 
In our model the 28/13 groundstate packing maximizes 
the number of coincident monomers from neighboring fil- 
aments, allowing a particular large number of favorable 
crosslinks to form. There are six monomers in a repeat 
unit of 28 monomers, which are (0,0); (4,1); (9,-1); 
(14,3); (19,1); and (24,-1) (see Fig. He)). 

Given a distribution of crosslinkers, we integrate out 
the spin degrees of freedom via a thermodynamic per- 
turbation theory. To lowest order this yields an effective 
Hamiltonian in terms of crosslinks alone, 



H 



eff 



l^2n t , m V(£,e!)n e , m >, (2) 



where e' = e — U/2. V {£,(.') is a pairwise coupling 
between cross-linking of different monomers along a fila- 
ment, 

V(£, £') = ^L. cos [wa(£ - £') - 2^(m - m')/6] e -l^'l/«* 

(3) 

Here £t = 2/3C is the twist persistence length, over which 
the oricntational correlations of the native filament ge- 
ometry are "washed out" by torsional fluctuations. This 
van der Waals-like coupling of distinct bonds reflects sta- 
tistical correlations in crosslinking along a filament. The 
rigidity of a crosslinking bond at layer £ pins the fila- 
ment in an orientation where certain nearby monomers 
are close to their most favorable binding direction, so that 
V {£,£') > 0. Hence, the range and strength of V (£,£') 
are determined by £t and U, respectively. 

The form of the effective binding model suggests that 
the statistics crosslinker of binding falls into the Ising 
or Bragg- Williams universality class. We analyze the 
mean-field thermodynamics within the grand canonical 
ensemble at fixed chemical potential, p,, which regulates 
the cost of removing a cross-linking protein from solu- 
tion. Assuming a constant mean site occupancy for the 
sites of the groundstate are occupied with a probabil- 
ity (n£_ rn ) = p, the mean-field equation of state is de- 
termined by the solution to the self-consistency condi- 



tion, p = (l + z 



-l„-up\ 



Here, z = exp[/3(/i + e )] 



is the effective fugacity of crosslinks, proportional to the 
concentration of unbound linkers in solution, and u is a 
measure of the net cooperativity of crosslinking. Specif- 
ically, u = N- 1 (3 J2' e ^ e V (£,£') where the sum is car- 
ried out over the total Nt, possible sites in the 28/13 



4 



groundstatc along a single filament. While coopcrativity 
monotonically increases with linker stiffness, u oc (f3U) 2 , 
this parameter has a more complex dependence on tor- 
sional rigidity. For small £ t , high-temperature, cooper- 
ative binding only occurs over short distances, so that 
u ~ £t- At larger values of £t the incommensurability 
between the native 13/6 and 28/13 twist symmetries re- 
quires significant distortions of either the filaments or 
the bonds between them. The incommensurate effects 
at long range lead to a reduction of u at large £t and 
maximum value around £t r* 60. 

The predicted mean-field equation of state is shown 
in Fig. [3j For low coopcrativity, u < u c = 4, p is a 
continuously increasing function of z, as crosslinking at 
distinct sites occurs largely independently in this regime. 
Increasing linker stiffness, increases the correlations in 
binding events, as indicated by rise in maximum linker 
susceptibility, Xp = dp/dz, for larger u. At the critical 
point u = u c , this susceptibility diverges at p = 1/2, Xp ~ 
\z — z c |~ 2 / 3 , indicating a second order transition. For 
u > u c , in which the stiff linkers enhances the cooperative 
interactions, the transition becomes first order with a 
discontinuous jump in linker density that increases with 
u. Owing to the Ising symmetry of Tie//, this model 
possesses a phase diagram for fixed z and U reminiscent 
of a "liquid- vapor" transition, in which a line of first order 
transitions terminates at a critical end point (see Fig. [3]) . 
Note that the value of the critical point implies a critical 
stiffness of order k c rts hsT/a 2 . 

A second result of this analysis is one-to-one correspon- 
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FIG. 3: (a)The predicted dependence of the mean site occu- 
pancy p on linker fugacity. The inset shows the phase diagram 
for fixed z and U. (b) The correlation between overtwist mea- 
sured by Im[ln g] and p for given values of £t. 



dence between mean occupancy of linker sites in bundles 
and filament overtwist observed in our SAXS measure- 
ments. We analyze the following monomer-monomer cor- 
relation function, g(£ ) = (exp {i(<j> to+se - (f> (o - lj 6£)}}, 
where the factor in the exponential is the excess angle 
between a monomer at £q and the next monomer £q + 51 
in 28-monomer packing relative to the native 13/6 twist, 
that is, the mean overtwist between neighboring bonds. 
We calculate g(£o) perturbatively to 0(U 2 ) for the given 
groundstate and find overtwist, as measured by Im[lng], 
to be continuously increasing function of p for any given 
values of linker and filament stiffness (see Fig. [3](b)). In- 
deed, because neighboring pairs of occupied bonds exert a 
torque on the filament to align monomers to the ground- 
state symmetry, it can be shown that Im[lng] ~ p 2 in 
the U — ► limit. Hence, the continuous (discontinuous) 
increase in crosslinking binding with increasing linker fu- 
gacity, implies a simultaneous continuous (discontinuous) 
structural transition in terms of filament twist. 

Theoretical results here suggest that F-actin crosslink- 
ing in parallel bundles is acutely sensitive to crosslinkcr 
flexibility. Both predicted regimes are experimentally 
observed. The continuous dependence of actin filament 
overtwist on the concentration of fascin, suggests that the 
these linkers are too flexible to exhibit a critical binding 
transition. While the comparative insensitivity of over- 
twist on linker concentration in espin bundles suggest 
that this binding occurs as a highly cooperative transi- 
tion, in which the rigidity of linkers immediately drives 
the bundle into the fully saturated and overtwisted state. 
The difference between espin and fascin binding sug- 
gests fundamental differences in the mechanism of bun- 
dle formation (twist, diameter, rigidity), which correlates 
to the distinct physiological functions of the respective 
actin bundles. Hair cells require structurally identical 
actin-bundles in order to mediate reproducible mechano- 
chemical transduction. This may be facilitated by an 
actin+crosslinker system in which the same bundle struc- 
ture is induced for a range of espin-actin molar ratios. In 
contrast, fascin's function is to organize cytoskeletal bun- 
dles in filopodial protrusions under a diverse set of me- 
chanical conditions [3] , a task that may be facilitated by 
the broad range of binding states and a sensitive depen- 
dence on the fascin-actin ratio. This view is consistent 
with in vivo observations of filopodial bundles that are 
weakly bound by fascin and highly dynamic 21 1 . 
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